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Abstract 

The resonance Z)*j(2317) which is considered to be the I^Pq state composed of charm 
and strange quarks has been discovered recently. The measured mass, which is about 
160 MeV lower than the mass of the I^Pq state obtained from the potential model 
calculation by Godfrey and Isgur, was considered surprisingly low and attracted a lot 
of theoretical investigations. We calculate the mass shift of the I^Pq state by using 
the coupled channel effect. Our result shows that the coupled channel effect naturally 
explains the observed mass of D*j(2317). 
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1 Introduction 



The BaBar Collaboration [1^ recently discovered a narrow resonance in D'^tt^, which 
is known as the D*j(2317), and it was confirmed by the CLEO ,2J and the Belle 
Collaborations Its decay patterns suggest a quark-model 0"*" classification. The 
measured mass, 2317.4±0.9 MeV 4j which is 40.9±1.0 MeV below the threshold of 
D^K~^, was considered surprisingly low compared to the predictions of the potential 
model calculations. For example, the prediction of the I^Pq mass by Isgur and Godfrey 
|S] was 2.48 GeV, and that by Eichten and Di Pierro [Hj was 2.487 GeV, which are 
about 160 and 170 MeV higher than the measured mass of D*j(2317). 

There have been many theoretical investigations which aimed to explain the mea- 
sured low mass of D*j{2317). Cahn and Jackson [7j used effective spin-orbit and 
tensor forces resulting from a model based on the nonrelativistic reduction. Bardeen 
et al. (8j worked with effective Lagrangians based on the chiral symmetry in heavy- 
light meson systems. Barnes et al. [HI considered a mixing between two molecular 
states \D^K~^ > and \D'^K^ >. Browder et al. ilO.^ proposed a mixing between the 
qq and 4-quark states and assigned a linear combination with less mass as D*j(2317). 

In this paper we calculate the mass shift of the I^Pq state composed of c and 
s quarks by using the coupled channel effect which was pioneered by Eichten et al. 
[TTj . We consider the coupling of the I^Pq bound state, whose bare mass is 2.48 GeV 
given in 0, with the two-meson continuum sectors of DK, D*K, DK* and D*K*. 
In reality, the D*K and DK* do not contribute since their spin-dependent statis- 
tical coupling coefficients are zero. We find that the contribution to the transition 
amplitudes from the Coulombic part of the interaction potential between quarks is 
about 60 % of that from the linear part. This aspect was first realized by Zambe- 
takis |12[ IT^ , and it is contrary to the usual consideration that the Coulombic part 
is negligible in the contribution to the transition amplitudes fTJ E]- As a conse- 
quence, the magnitude of the reduced interaction hamiltonian resulting when both 
linear and Coulombic parts are included is about 2.6 times that resulting when only 
the linear part is considered. This fact is crucial in the analysis of this paper: When 
the Coulombic part is not included, there is no physical mass eigenstate below the 
threshold of D^K^ for physically acceptable parameterizations. On the other hand, 
when the Coulombic part is included, there exists a physical mass eigenstate with its 
mass eigenvalue close to the measured D*j(2317) mass in a natural manner. This 
paper is organized as follows. In section 2 we briefiy summarize the coupled channel 
effect. Even though it is explained in detail in Refs. El we think that a 
concise summary is useful for this paper. In section 3 we explain the procedures and 
results of our calculation. The last section is conclusion. 



2 Coupled Channel Effect 
2.1 Basics 

The coupled channel effect takes into account the fact that two-meson continuum 
states of DK couple to two-quark cs bound states. There are two important effects 
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of this couphng. The first is to give rise to mass shifts and configuration mixing of 
cs bound states. The second is to provide broad peaks of cs resonance states above 
threshold |TTl IT^ IT^ . In this section we briefly summarize the coupled channel effect 
which has been explained in detail in Refs. [TTJ [T21 HH] • 

When we take the coupled channel effect into consideration, Hilbert space is com- 
posed of two sectors: discrete cs bound states and continuum states of two mesons. 
Then the total hamiltonian is written as 

H = H, + H,c, where H,^("^^ ^ J , ^pc = f ( ^° _^ Y)' « 

Hq and He are diagonal respectively in the quark- antiquark bound state space \i > 
and in the two-meson continuum space |P,A > (where P is the relative momentum 
of two mesons), which satify 

i/oK >= M^'-^lz >, Hc\P,\>=EciP)\P,\> , (2) 

<i\j>=Sij, <P,A|P',A'>=(27r)353(P-P')5AA', <i|P,A>=0. (3) 

Hpc in (0), where PC means pair creation, connects two spaces \i > and |P, A >. 

Physical mass eigenstates are eigenstates of the total hamiltonian if in (P). We 
denote the physical mass eigenstate by |iV > with a collective quantum number N. 
Then with a physical mass M^r, \N > satisfies 

H\N>=Mn\N>. (4) 

When the value of ^ in (P) varies to zero, > becomes one of the bare states \i >. 
Physical mass eigenstate > can be expanded in terms of bare states as 

|iV>=^a,|z>+^ J J^sh{P)\P,X>. (5) 

From ((I}-© bx{P) and in © are related by 

fp\ V- < P,MHpc\i > .^x 
^'^'^^=^ M,-Ec{P) 

then from and (jlj) we get 

H,sa^ = MNa^ with H^s = M^""'" + n{MN), (7) 

where M^^^^ is a diagonal matrix with M}'^^^ in (j2I) as diagonal elements, is a 
column vector with probability amplitudes af as its elements, and the matrix elements 
of fl{W) is given by 

O (W) V / <^|gpc|P,A><P,A|gpc;|j> 
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For W above threshold, n{W) has the structure of n{W) = A{W) - iT{W)/2 with 
A(H^) and r(iy) real symmetric matrices. Once the matrix fi(iy) is known, the mass 
eigenvalues of the coupled system are found from 

Det(M^J - [M^^'^ + Q{Mn)]) = 0. (9) 

We can calculate the normalization constant of the amplitudes in (jS)) from 

A/''(El«^r + E/|^l&A(P)P) = i. (10) 

Using (jH)) and ((HI), we get 

Af' = {J2\a,\^ + J2<^iAMN)^,y\ (11) 

where 

u;,,{M^) = -^Q^^{W)\w=M„. (12) 

The probability that the physical mass eigenstate > be in the quark- antiquark 
bound state sector is given by Z = N"^Y^i Idjp, while that in the two-meson state 
sector is 1 — Z. 



2.2 Mechanism of Light Quark Pair Creation 

The Cornell group studied the effect of OZI allowed decay channels ^1]. They pro- 
posed that the following interaction hamiltonian is responsible for the decay as well 
as the binding of quark-antiquark bound states. 

Hi = J d'^d'y ■ Pa(x)\/(x - y)p,(y) :, (13) 



a=l 

where 

'[/(^\ - _1 J 

-.2- 



y(r) = -- + 4. (14) 



In ffT^ Pa(x) = ip"^ {'x)^Xaip{'x.) are the color densities of quark fields. This model 
corresponds to the vector coupling since (jl3|) is the leading term of the vector coupling 
hamiltonian in the nonrelativistic expansion. 

One can calculate the transition amplitude between two-quark bound and two- 
meson continuum states using the interaction hamiltonian (|13p. Then (jS)) becomes 



nL,mL' 



(W) = Y.j P^dP 



W - E,{P) - E2{P) + te 



= T.lKL,mL'm - l^L^n^AW)], (15) 



4 









D*+ 


K+ 






1864.6 


1869.4 


2006.7 


2010.0 


493.677 


497.648 


891.66 896.10 



Table 1: Meson Masses (MeV) used in the calculation jlj. 



where 

HL,mL'iP) = fY.C{Js,LL\JM)lUP)lLL'{P) (16) 
I 

with 

2 1 

f = = — • (17) 

In fll6|) . lli{P) is the momentum dependent factor of the transition amplitude. The 
superscript / in (fT3j) denotes the participating coupled channels which are composed 
of two mesons. The angular momenta and energies of these two mesons are denoted by 
Ji, J2, and Ei{P), E2{P), respectively. H^^^ ^i^,{P) depends on the coupled channel, 
and it also depends on J and s (the total angular momentum and spin of the quark- 
antiquark bound state). 

For the vector coupling, /^^(P) in (fT?)|l is given by 



with 



e(t) = [te-^' + {t^-l)e-^'/'^erf{^)]+4Pa\[-te-^\e-'^ (19) 

where Rniir) is the radial wave function of the quark- ant iquark bound state, jiit) is 
the spherical Bessel function, and pq = mQ/{mq + mq). The first and second terms 
in ()19|) come from the linear and Coulombic parts of (fT^ . respectively. 

3 Mass Shift of D*j (2317) Meson from Coupled 
Channel Effect 

Isgur and Godfrey obtained the mass of the l^Po state composed of c and s quarks as 
2.48 GeV from the interaction between two quarks. Their result corresponds to the 
bare mass M^^^'^ in Q. In this section we calculate the mass shift of the I^Pq state 
using the coupled channel effect explained in the previous section. In particular, ^ 
gives the physical mass given when the coupled channel effect is incorporated. 
For V{r) in ()14|1 we use the following two parameterizations in our calculation. 

(A) Eichten et al. |llj: k=0.517, a=2.12 GeV, mc=1.84 GeV. 

(B) Hagiwara et al. k=0A7, a = 1/V(U9 GeV, m^=1.32 GeV. 

For the masses of s and u, d quarks, we use ms=0.55 GeV, mu^rf=0.33 GeV. 
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P (GeV) 



Figure 1: /n(P) for the potential A with /ib=0.46 GeV. 



FoUowing we use the radial wave functions of cu, cd, su, sd states with L = 

3 1 

which are given by the ground state harmonic-oscillator wave functions (2/1^2/714) 
exp(— and that of cs state with L = 1 by the first excited harmonic-oscillator 
wave function (Sa/i^s/SsTri) r exp(— /if,^'^^/2). P in (fT7|) and /i^ here are related by 
(3 = fJ^af^- the value of (3 in the calculations of (fTTjl and (fTH|) . we take /ia=0.61 GeV 
which is the average value of fia{K)=0.61 GeV and /ia(-D)=0.60 GeV given in [TB] . 
For the value of pg in (fT8|l defined as mQ/{mq + mg), we take mq as (mc + ms)/2. 
The formula for /^^(P) in jlHI) was derived for the case where two mesons in the 
coupled channel have the same heavy quark such as DD [llj. However, in our present 
study two mesons are D and K mesons, and ()18j) is not exact in our system. So, we 
approximate our system by using the formula for P^j^{P) given in ()18j) with the above- 
mentioned values of f3 and pq. For the value ph which comes into the cs state radial 
wave function Rnbij) in ((IHI), we take two prototype values: /ife=0.46 GeV given 
in jin] for i5sj(2632) and /ib=0.57 GeV in [17J for ^Po(cs). In the calculation in 
this paper we also adopt the approximation of taking only one cs bound state of 
I'^Po for the quark-antiquark bound state sector, whereas we take all the two-meson 
continuum sector of DK, D*K, DK* and D*K*. That is, we do not include in our 
calculation the higher cs bound states of h^Pq whose masses are located above the 
threshold. If we include those higher cs bound states in the calculation, the physical 
mass eigenvalue would be modified. However, we expect the modification is not large 
since the I^Pq is very close to the threshold and then the coupling channel effect to 
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Figure 2: Ifi{P) for the potential A with /ib=0.46 GeV. 

this state is dominant. There is also a merit of this approximation that one can see 
clearly how the coupled channel effect gives rise to the mass shift. 

The statistical coupling coefficient C( Js, LL', Ji J2, /) in (fT^ is 1 for / = in 
the DK channel, 1/3 for / = and 8/3 for / = 2 in the D*K* channel, and other 
coefficients are zero. In the calculation we use the meson masses given in Table 1 |3]. 
We calculate /^^(P) using (fTHj) . The results of /^^^(P) for n = L = 1 and / = 0, 2 are 
presented in Figs. ^andlH As we see in the figures, the contribution to /^^(P) from 
the Coulombic part is about 60 % of that from the linear part. 

Ref. ^] ignored the Coulombic part of V{r) in (fT^ in order to simplify their cal- 
culations. They argued that this is justified because small quark separations are not 
important in hadronic decays. When Eichten et al. jH] analyzed the recently found 
X(3872) state using the coupled channel effect, they also did not include the Coulom- 
bic part in (fT!^ in their calculation by considering that it is a good approximation 
to neglect all effects of the Coulomb piece in (fTn|) . It has been a common practice to 
ignore the Coulombic part in (fT!^ . However, it was first found by Zambetakis [T^IT!^ 
that if one includes the Coulombic part in (fT!?|) . he finds a significant modification of 
ll^i^. In fact, we find in this paper that the contribution to P^j^ from the Coulombic 
part is about 60 % of that from the linear part, and then the magnitude of VL resulting 
when both linear and Coulombic parts are included is about 2.6 times that resulting 
when only the linear part is considered. 

The fact that the inclusion of the Coulombic part in (|T^ increases the magnitude 
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Figure 3: The result when the potential B and fib= 0.46 GeV were used. The upper 
first and second lines are fl when the Coulombic part is not included and when it 
is include, respectively, and the lower diagonal line is W — M^''^'^'^ where M^^^'^=2A8 
GeV given in 0. As we see here, fl without the Coulombic part does not meet 
with W — M^^'^'^, which means that there does not exist a solution of the eigenvalue 
equation Q below thresholds. On the other hand, Q with the Coulombic part meets 
with W — M^'^^^ below thresholds and the value of 14^=2.31 GeV at the meeting point 
is the eigenvalue Mjy of ©• 
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Parameterizations 


Mn 


AMtv 


Z 


Potential 


^ib (GeV) 


(GeV) 


(GeV) 




A 


0.46 


2.29 


-0.19 


0.67 




0.57 


2.27 


-0.21 


0.70 


B 


0.46 


2.31 


-0.17 


0.67 




0.57 


2.29 


-0.19 


0.71 



Table 2: The results of M^r, AMi\i and Z for each parameterization of the potential 
and the fib value, when Pq(= mQ/{mQ + rriu^)) with mq = {nic + rn^jl was used. 



Parameterizations 


Mn 




Z 


PQ 


Potential 


Pb (GeV) 


(GeV) 


(GeV) 






A 


0.46 


2.31 


-0.17 


0.63 


Pc 




0.57 


2.29 


-0.19 


0.68 




B 


0.46 


2.32 


-0.16 


0.62 






0.57 


2.31 


-0.17 


0.67 




A 


0.46 


2.23 


-0.25 


0.73 


ps 




0.57 


2.19 


-0.29 


0.75 




B 


0.46 


2.27 


-0.21 


0.73 






0.57 


2.25 


-0.23 


0.75 



Table 3: The results of M^, AM^r and Z for each parameterization of the potential 
and the /ib value, when Pq{= mQ/{mQ + niu^)) with mg = rric or was used in 
order to see the sensitivity of the results to the parameter values. 



of Vt about 2.6 times is very important in the analysis of this paper. Fig. |21 is the 
result which we get when the potential B and /ib= 0.46 GeV are used. The upper 
first and second lines are VL when the Coulombic part is not included (that is, when 
only the linear part is considered) and when it is included, respectively, and the lower 
diagonal line isW — M^^^'^ where M^'^'''^=2.48 GeV given in [S]. As we see the figure, 
VL without the Coulombic part does not meet with W — M^^'^'^ , which means that 
there does not exist a solution of the eigenvalue equation © below thresholds. On 
the other hand, VL with the Coulombic part meets with W — M^^^'^ below thresholds 
and the value of 1^=2.31 GeV at the meeting point is the eigenvalue of 
The experimentally measured value of D*j(2317) mass is 2317. 4±0. 9 MeV, which is 
40.9±1.0 MeV below the lowest threshold energy 2358±0.5 MeV of D^K+. For all 
the cases of the parameterizations given in Table 2, we have similar situations as that 
shown in Fig. El In Table 2 we present the physical mass eigenvalue and the 
mass shift AM^v = Mm - M'^^=^*^(2.48 GeV). The values of Z which we calculated 
using (fTTj) are also presented in the table. 

In order to see the sensitivity of the results to the value of the parameter pq, 
we also performed the same calculation for pq = p^ and p^, where mq = rric and 
vis in Pq = mq/{mq + rriu^d), and the results are presented in Table 3. For all 



9 



these different parameterizations, it is still true that the contribution to /^^^ from the 
Coulombic part is about 60 % of that from the linear part, and then the magnitude 
of Q resulting when both linear and Coulombic parts are included is about 2.6 times 
that resulting when only the linear part is considered. Among all the cases in Table 2 
and 3, only for the cases in the 5th and 6th rows of Table 3 there exists an eigenstate 
also for the Q with k = 0, with its eigenvalue which is just belw the threshold of 
D^K^ . However, the parameterizations of those cases are not physically acceptable 
since their mass eigenvalues presented in the 5th and 6th rows of Table 3 are too 
small compared to the measured value. 

4 Conclusion 

It has been a common practice to ignore the Coulombic part in the calculation of 
the transition amplitude between two-quark bound and two-meson continuum states. 
However, the calculation in this paper shows explicitly that it is about 60 % larger 
when both the linear and Coulombic parts are included in the calculation compared to 
the case where the latter part is ignored. Consequently, the magnitude of becomes 
about 2.6 times larger when the Coulombic part is also included. This fact is crucial 
in the analysis of this paper: When the Coulombic part is not included, there is 
no physical mass eigenstate below the threshold of D'^K^ for physically acceptable 
parameterizations. On the other hand, when the Coulombic part is included, there 
exists a physical mass eigenstate with its mass close to the measured L'*j(2317) mass 
in a natural manner. Therefore, it was shown that the coupled channel effect explains 
why the mass of D*j(2317) is about 160 to 170 MeV lower than the bare mass of the 
l^Po bound state obtained from the most potential model calculations. 
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